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Extended Self-Similarity (ESS), a procedure that remarkably extends the range of scaling 
for structure functions in Navier-Stokes turbulence and thus allows improved determi- 
nation of intermittency exponents, has never been fully explained. We show that ESS 
applies to Burgers turbulence at high Reynolds numbers and we give the theoretical 
explanation of the numerically observed improved scaling at both the infrared and ul- 
traviolet end, in total a gain of about three quarters of a decade: there is a reduction of 
subdominant contributions to scaling when going from the standard structure function 
representation to the ESS representation. We conjecture that a similar situation holds for 
three-dimensional incompressible turbulence and suggest ways of capturing subdominant 
contributions to scaling. 

o 

CN ' 1. Introduction 



Extended Self-Similarity (ESS), discovered bv lBenzi et al. (1993), is the empirical ob- 



ov 

servation that in fully developed turbulence, when plotting structure functions of order 
p vs, say, the structure function of order three, rather than the traditional way where 
they are plotted vs the separation, then the range over which clean power-law scaling is 
' observed can be substantially increased. This has allowed a much better determination 

of the scaling exponents Cp of the structure functions of order p — or at least of ratios of 
such exponents — and has been key to confirming that three-dimens i onal h igh Reynolds 



number incompressible turbulence does not follow the iKolmogorovl ([19411 ) scaling laws 
Cp = J*/3, but instead has anomalous scaling, whose exponents cannot be obtained solely 
through dimensional arguments. 

In spite of several attempts to explain the success of ESS (see, e.g.lBhattachariee fc SairJ 
1999; Sa in fc Bhattachariedll999tlFuiisaka fc GrossmannlbOOluSegel et al.lll996tlYakhotl 
200l[ and Section [5]) , the latter is still not fully understood and we do not know how 



much we can trust scaling exponents derived by ESS. It would be nice to have at least 
one instance for which ESS not only works, but does so for reasons we can rationally 
understand. A very natural candidate might be the one-dimensional Burgers equation. 
Early attempts to test ESS on the Burgers equation did not show any appreciable in- 
crease in the quality of scaling through the use of ESS. As we sh all see in Section EE the 



conclusion that "ESS does not work for the Burgers equation" (jBenzi et al.l 119951 ) was 
just reflecting the computational limitations of the early nineties. 

In Section [2] we recall some basic facts and notation for ESS in three-dimensional 
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Navier-Stokes turbulence. Then in Scction[3]we turn to the Burgers equation and present 
new numerical evidence that ESS works for Burgers, provided high enough spatial resolu- 
tion is used. In Section[3]we use asymptotic theory to explain in detail why ESS works for 
the Burgers case. Finally, in Section [5] we examine the possible lessons from our Burgers 
ESS study for three-dimensional Navier-Stokes turbulence. 



2. ESS in a nutshell 

Consider the three-dimensional Navier-Stokes (3DNS) equation 

d t v + v ■ Vv = -Vp + ^V 2 v, V-v = 0. (2.1) 

For the case of homogeneous isotropic turbulence, (longitudinal) structure functions of 
integer order p are defined as 

S p (r) = ({Sv l{ (r)) p ), 
in terms of the longitudinal velocity increments 

r 

r 



6vu(r) = [v(x + r) — v(x) 



(2.2) 



(2.3) 



where r = |r| and the angular brackets denote averaging. There is experimental and 
numerical evi dence that, at high Reynolds numb ers, structure functions follow scaling 
laws (see. e.g. iMonin fc YariomllT97llFrischlll995l) 



S p (r) oc r^ p 



(2.4) 



over some range of separations (the inertial range) L ^> r 3> rj p . Here L is the inte- 
gral scale and rj v the dis s ipation scale. The latter m ay depend on the order p (see e.g. 
Paladin fc Vulpiani|[l987t iFrisch fc Vergassolalll99llh 

Of course, the dominant-order behaviour given by (|2.4[) is accompanied by subdom- 
inant corrections involving the two small parameters characteristic of inertial-range in- 
termediate asymptotics, namely r/L and T) p /r. The simplest would be to have 

S p (r) = C p r*> (l + Df{r/L)C + D^Wr)*™) + h.o.t. , (2.5) 

where h.o.t. stands for "higher-order terms" and where g p R > and g p v > are the 
infrared (IR) and ultraviolet (UV) gaps, respectively. For a given Reynolds number and 
thus a given ratio L/rj p , the smaller the gaps and the constants D p R and D p v , the larger 
the range of separations over which subdominant corrections remain small. 

The ESS is an operational procedure that effectively enlarges the range of separations 
over which dominant-order scaling is a good approximation. In its simplest formulation, 
one considers two integer orders n and m and plots |5„(r)| vs |S m (r)| and finds empirically 
that the scaling relations 

\S n (r)\^\S m (r)\ a ^ m \ (2.6) 

with suitable exponents a(n,m), hold much better than (|2.4[) . One particularly interest- 
ing instance of this procedure is when m — 3. We then kno w from Kolmogorovl (|1941r ) 
that, to dominant order, we have the four-fifths law (see, also Frisch 1995h 

S 3 (r) = ~er, (2.7) 



where e is the mean energy dissipation per unit mass. Thus, the third-order structure 
function (divided by — (4/5)e) may be viewed as a deputy of the separation r. A variant 
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of the ESS, which frequently gives even better scaling, is to use alternative structure 
functions, defined with the absolute values of the longitudinal velocity increments, namely 

F p (r) = <|^||(r)| p ) . (2.8) 

It is then found empirically that 

F n (r) w F m (r)^"' m > , (2.9) 

with suitable scaling exponents /3{n, m). Whatever its empirical merits, the variant proce- 
dure has the drawback that there is no equivalent to the four-fifths law for the third-order 
structure function with the absolute value of the longitudinal velocity increment. Thus 
we cannot safely use F 3 (r) as a deputy of r. We shall come back to this in Section [5] 

3. ESS revisited for the Burgers equation 

The one-dimensional Burgers equation 

dtU + ud x u — vd^,u; u(x, 0) = uq(x), (3-1) 

whi ch was introdu ced originally as a kind of poor man's Navier-Stokes equation (see, 
e.g. . lBurgersl 1974 ). has some dramatic differences with three-dimensional N avier-Stokes 



(3DNS) turbulence, foremost that it is integrable (|HopJll950t ICold 119511) and — as a 



consequence — does not display self-generated chaotic behaviour. Nevertheless i t does 



display anomalous scaling in the following sense: superficially, the iKolmogorovi (|194lh 
theory is applicable to the Burgers equation as much as it is to 3DNS. However, when 
starting with smooth initial data uq(x), the evolved solution in the limit of vanishing 
viscosity v will display shocks. Thus ( p = 1 for p > 1. When the Reynolds number is 
finite, structure functions will display scaling only over a limited range of separations r. 
Therefore, the Burgers equation may be a good testing ground for ESS and also perhaps 
for understanding why and when it works. Such considerations did not escape the creators 
of the ESS technique. Unfortunately, no clean scaling for structure functions is observed 
with the Burgers equation, either in the standard representation or in ESS, as long as 
simulations are done with the spatial resolution easily available in the early nineties, 
namely a few thousand collocation points. Scaling emerges only at much higher spatial 
resolutions with 128K (128 x 1024) Fourier modes and becomes fully manifest with 256K 
modes, which is now also the highest resolution achievable numerically within a time 
span of a few week. 

Let us now explain our numerical strategy for studying ESS with the Burgers equation. 
Our goal in doing preliminary numerical experiments is to understand ESS in a rational 
way, starting from the basic equations. For this it is advisable to keep the formulation 
minimally complex, avoiding too many "realistic trappings" . For example, there is no 
need at first to assume random initial conditions: we can just take an L-periodic initial 
condition and define the structure functions by integrating over the period: 

S p (r) = (1/L) [ dx[u(x + r,t) -u(x,t)] p . (3.2) 



We shall mostly work with a very simple single-mode model for which the initial condition 
is 27r-periodic deterministic and has a single Fourier mode 

ito = sinx. (3-3) 

As we shall see in Section 01 it is easy to extend the theory from the deterministic to 
the random case. 




Figure 1. Compensated sixth-order structure function in standard (continuous line) and ESS 

(dashed line) representations. 



We integrated the Burgers equation (|3.1|) with the initial condition (|3 . 3|) using a 
pseudo-spectral method with N collocation points and a two-thirds alias-removal rule. 
Time stepping was done in double precision by a 4th order Runge-Kutta scheme with 
constant time step St. The viscous term was ha ndled by the slaving technique known as 
ETDRK4 described in ICox fc Matthews! (|2002h . which allows taking a time step about 



ten times larger th an would be permitted with a direct handling of the viscous term. (It 
was pointed out bv lKassam k. Trefethen ( 2005 ) that ETDRK4 can produce numerically 



ill-conditioned cancellations; following a suggestion of J.Z. Zhu (private communication), 
we handled these by performing Tayl or expansions to suitable order rather then by the 
complex-plane method proposed bv lKassam fc Trefethenl ( 20051 ).) Slaving results not only 



in considerable speed-up but in much less accumulation of rounding noise. 

The parameters of the run were N = 256K and St = 10~ 5 . Output was processed at 
t = 2, well beyond the time of appearance of the first shock at t = t* = 1 . 

Fig. [T] shows the compensated structure function — that is divided by the theoretically 
predicted inertial-range dominant term — of order six in both the standard representation 
and in a variant of the ESS representation. Our variant uses 

Ss(r) = (3.4) 

where e = —(\/L)(d/dt) J Q L dxu 2 /2 is the mean energy dissipation. It is e asy to show that 
the B urgers counterpart of the four-fifth low is a "minus twelve" law ( cf. iGurbatov et aT 



1997T ) which makes S3 the appropriate deputy of the separation r. In our opinion it is 
important to chose the constant in the definition of S3 in such a way that it becomes r 
with a unit factor (to dominant order). Otherwise an ESS plot in log-log coordinates may 
show an overall improvement in quality of scaling, without our being able to disentangle 
the small-separation (UV) improvement from the large-separation (IR) improvement. As 
we shall see, both are in general present and have quite different origins. 

Fig. [T] shows a substantial improvement in scaling for ESS, that is a wider horizontal 
plateau in the compensated structure function; and this at both the IR and UV ends. 
This is the first evidence that "ESS works for Burgers". Next we shall understand why 
it works. 
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4. Asymptotic theory of ESS for the Burgers equation 

We now give the theory for improved ESS scaling when p > 3, first for the single mode 
case and then upgrade it for the case of random solutions. 

To handle the infrared (IR) contributions to the structure functions we can work with 
an infinitely sharp shock, taking v — > 0. The dominant contribution to structure functions 
of integer order comes clearly from intervals [x, x + r] which straddle the shock location 
x s (in this Section the time variable is written explicitly only when needed). It is also 
easily shown that for p > 3 the first-order subdominant contributions comes from the 
small changes of the velocity, immediately to the left and the right of the shock, which 
are expressible by Taylor expan ding the velocity to first order in these two regions (see 
Section 4.2.2 of lBec et al.ll2000h : 

u(x) = u_ + (x — x s )s- + h.o.t., u(x) — u + + (x — x s )s + + h.o.t., (4-1) 

where u- and u + are the velocities immediately to the left and to the right of the shock 
and s_ and s+ their respective gradients. Starting from (13. 2|) and limiting the integration 
domain to the interval [x s — r, x 3 ], which corresponds to the straddling condition, we 
obtain, using (|4.ip 



LS p (r) = {-l) p (a*Y - A^ 1 ^ + s-)|r 2 ) + h.o.t., (4.2) 

where 

A = u- - u+ > (4.3) 

is the amplitude of the shock and L = 2ir is the spatial period. Specialising to the 

third-order structure function and to its rescaled version the separation deputy S 3 (r) = 
Sz(r) / (—12e) , we obtain 

LS 3 {r) = -A 3 r + Q \ A 2 (s_ + s+)r 2 + h.o.t. (4.4) 

&(r) = r - {a- + s + )r 2 + h.o.t., (4.5) 



2A 



where we have used the relation 



A 3 

e = —— (4.6) 



12L 

between the energy dissipation and the shock strength. 

We now eliminate r between (|4.5I) and (14. 2|l . so a to rewrite the structure function of 
order p as an expansion in the separation deputy : 

LS P = (-l) p (A p S 3 - AP- 1 is- + s+)Si) + h.o.t. (4.7) 



Comparison of the "standard" expansion (|4. 2|) of the structure function and its ESS 
expansion ()4.7|) . shows that they have the same dominant terms and that their first 
subdominant corrections differ only by a numerical coefficient: p/2 for the standard case 
and (p — 3)/2 for ESS. Hence the subdominant correction has been decreased by a factor 
p/(p — 3). For the case of the sixth-order structure function, considered in Fig. [I] this 
is a reduction by a factor two. Hence, the ESS inertial-range scaling extends by a factor 
2 further into the IR direction before a same level of degradation is achieved as for the 
standard case. Note that for large ps the gain in scaling range becomes smaller. 

Next we turn to the ultraviolet (UV) contributions which now require a finite viscosity 
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v that broadens the shock. Standard boundary layer analysis for the shock in the frame 
where the shock is at rest (which basically amounts to dropping the time derivative term 
in the Burgers equation) gives the following well-known tanh structure : 

u(z)«--tanh— . (4.8) 



2 4i/ 

Hence the UV structure functions are given by 



xA t (x + r)A 



tanh tanh 



Av Av 



(4.9) 



We are interested in the expansion of these structure functions for r much larger than the 
typical width Av/ A of the shock. For this, we first established the following expansion 
(for large f) 

/>oo 

dx [tanh(x) - tanh(z + f)] p = (-2) p (f - H p _x) + t.s.t. (4.10) 

i 

where "t.s.t." stands for "transcendentally small terms" such as exp(— f) and 

H p = 1 + 1/2 + 1/3 + ... + 1/p is the Harmonic function, which behaves as lnp for large 

p. Using (|4.1ip in l|4.9p . we obtain 

LS^r) = (-l) p (A p r - 4i/ J ff p _ 1 A p - 1 ) + t.s.t. (4.11) 

Specialising to p = 3, we have 

LS 3 = -A 3 r + 6vA 2 + t.s.t., S 3 (r) = r - 6-^ + t.s.t.. (4.12) 

Proceeding as in the IR case, we re-expand S p in terms of the deputy separation S3: 

LS p (r) - (§ 3 AP - Av (h p -i - A^" 1 ^ + t.s.t. (4.13) 

Thus, we see that with ESS the subdominant UV term for the structure function of order 
p is reduced by a factor 2i7 p _i/(2iJ p _i — 3) and, again, the range of scaling is extended 
by the same factor. For p = 6 the extension factor is 137/47 2.91. Combining the UV 
and the IR gains, we see that the scaling for Sq is extended by a factor 5.92, that is about 
three quarters of a decade. 

Next we upgrade the arguments to the case of the Burgers equation with smooth 
random initial conditions and forcing, defined on the whole real line. We assume that 
the forces and the initial conditions are (statistical) homogeneous and have rapidly de- 
creasing spatial correlations (mixing) . We can then use ergodicity to obtain the following 
representation of structure functions: 

1 ' L 



S p {r) ee ((u(x + r) ~u{x)) p ) = lim - / dx (u{x + r) - u(x)) p , (AAA) 

L^oc L J Q 

where the limit is in the almost sure sense. In the present context we have typically an 
infinite number of shocks on the whole line, but a finite number per unit length. The 
shock amplitudes A and the left and right velocity gradients s_ and s + become random 
variables. Revisiting the arguments given above for the deterministic single-mode case, 
we find that in both the IR and UV expansions we now have to add the contributions 
stemming from the various shocks. Using ergodicity we obtain, in the UV domain 

S p (r) = «A*> r - AvHp^iAP- 1 )) + t.s.t. (4.15) 
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and, in the IR domain 

S p (r) = (-If ((A p ) r - (Af- 1 (s+ + |r 2 ) + h.o.t. (4.16) 

In the UV case we now make use of the inequality (A P )/(A P_1 ) > (A 3 )/(A 2 ), which 
follows, for p > 3, from the log-convexity of the moment function q ln(A 9 ) for a positive 
random variable A. We can then finish the analysis of the depletion of subdominant 
contributions essentially as done above (in the UV) case and conclude that ESS extends 
the UV scaling by at least a factor 2H p _i/ (2H p _i —3). In the IR domain, the presence of 
the random quantity s + + s_ prevents us from reaching similar conclusions and it is thus 
not clear that there is a general result regarding improved ESS scaling in the IR regime. 
This can however be circumvented, if we assume (i) that the Burgers turbulence is freely 
decaying (no forcing) and (ii) we limit ourselves to times large compared to the typical 
turnover time of the initial condition. The solution degenerates then into a set of ramps 
of slope exactly 1/t, separated by shocks. Hence s+ + s_ « 2/t, which is deterministic. 
Thus, using the same log-convexity inequality as above we infer that ESS extends the IR 
scaling by at least a factor p/(p — 3). 



5. Back to three-dimensional Navier Stokes turbulence 

Here, we have explained the success of ESS by a depletion of subdominant IR and UV 
contributions. How does this relate to various explanations given over the past one and 
half decade for 3DNS turbulence? Let us mention a few . Sain and Bhattacharjee (see, 
Bhattachariee fc Sain 19991 : Sain fc Bhattacharieei 19991) resorted to phenomenology in 



proposing cross-over functions from the dissipation- to the inertial-range for the structure 
functions defined in Fourier space. They inferred that in the UV regime, the linear scaling 
in log-log plots of structure f unction deteriorate at lower v alues of wavenumbers than in 
the corresponding ESS plots. iFujisaka fc Grossmannl (|2001l ) — again phenomenologically 
— introduced a scaling variable to include crossovers between various subranges of scaling 
behaviour for the magnitude of longitudinal velocity differences and thereby showed that 
the scaling improves at the IR end when ESS is used. In the spirit of our formulation of 
a theory for ESS, the work of Segel et al. ( 19961 ) comes closest. In this paper they have 
correctly identified the mechanism of increased scaling range at the UV end: a reduced 
coefficient in the subdominant diffusive correction. Note that their work on ESS was in 
the context of inte rmittency and an o malou s scaling for passive scalar dynamics. They 
used the mo del of Kraichnan (| 19681 . 119941) and supplemented it by a closure relation 
suggested by iKraichnan ( 1994|). which was later shown not to be fully consistent with 



the original model (see e.g.. iFalkovich et al.ll200lh . Benzi (private communication) also 



made a similar observation using a passive scalar shell model. 

How much of our findings for the Burgers equation carry over to 3DNS? It is important 
to realize that the improved scaling can be both at the IR and the UV end of the scaling 
regime. In order to avoid mixing up the two types of improvements one should correctly 
calibrate the choice of the deputy for the separation r, by using the four-fifths law for 
3DNS: 

If a given turbulent flow shows a gain in scaling when using ESS, either in the IR or the 
UV domain or both, we can then use (|2.5[) to interpret the gain in terms of subdominant 
corrections. It is best to handle the IR and UV cases separately. Following the same 
procedure as in Section [4] it is easily shown that ESS gives just a modification of the 
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coefficient of the first subdominant contribution provided the gap g p between dominant 
and first subdominant exponent is independent of the order p. The fact that ESS works so 
nicely for 3DNS suggests that this independence may actually hold. If so, it is immediately 
seen that the reduction in the subdominant coefficient is equal to the gain in scaling raised 
to the gap value. 

We begin to see here ESS as a way to obtain information on subdominant corrections. 
This is of interest for several reasons. For e xample, subdominant corrections ca n give 
rise to spurious multifractal scaling (see e.g. jAurell et al.lll992t iMitra et al.ll2005l ). Fur- 
thermore, consideration of subdominant corrections is needed to expla in the absence of 
logarithms in the third-order structure function (cf. iFrisch et al.ll2005l) . Also, the multi- 
fractal description of turbulence is quite heuristic and arbitrary and would be much more 
strongly constrained if we had information on subdominant terms and on gap values. 

This may be the right place to discuss the issue of the best deputy of separation in the 
ESS procedure. Should one use Ss(r) or the function Fs(r) that is defined with the third 
moment of the absolute value of the longitudinal velocity increment and which is easier 
to extract from experimental data because it involves only positive contributions? It is 
not clear to what extent the two procedures are equivalent. For the case of the Burgers 
turbulence it is easy to show that — 5a(r) and F 3 (r) have the same dominant and first- 
order subdominant terms: they differ only in subsubdominant contributions. For 3DNS 
longitudinal velocity increments are somewhat more likely to be negative than positive 
but they have no reason to have exactly the same scaling behaviour alth ough the scaling 
exponents are found to be nearly equal (jVainshtein fc Sreenivasanlll994l ). There has also 
been some amount of discu ssion regarding th e scaling of longitudinal and transverse 
structure functions (see e.g.. iBenzi et al. 20091 ). It may well be that they differ only by 
the relative strength of subdominant terms. 

We finally address the issue of appropriate strategies to systematically obtain sub- 
dominant corrections from experimental or simulation data. The most straightforward 
method is to determine the dominant-order contribution using ESS and then to sub- 
tract it from the the data. The result can then be analysed either in the standard way 
or in the ESS fashion. However it is known that when subdominant terms are rather 
sizeable it is better to determine them at the same time as the dominant ones. One 
instance is the simultaneous determination of dominant-order isotropic scali ng and sub- 



dominant anisotropic corrections for weakly anisotropic turbulence (see e.g.. lArad et al 
19981 : Biferale fc Procaccia 20051 and references therein). 



The determination of both dominant and subdominant terms can be much improved 
when higher precision is available, such as is frequently the case in double-pre cision 
spectral simulations. Indeed it was recently pointed out by Ivan der Hoeven (2009) that, 
when trying to extract the asymptotic expansion of a function f{x) as x — > oo from data 
sampled at a large number of values of discrete x- values, it not advisable to first try to 
obtain the dominant order and only then to look for subdominant terms: without the 
knowledge of subdominant corrections, the parameters appearing in the dominant term 
will be very poorly conditioned. It is better to completely subvert the dominant-order- 
first strategy by introducing the method of asymptotic interpolation, which relies on a 
series of transformations, recursively, peeling off the dominant and subdominant terms in 
the asymptotic expansion without any need to know their detailed functional form. These 
transformations can, in principle, be carried out till either the rounding noise becomes 
significant or asymptoticity is lost. At the end of the process, the newly transformed data 
admit a simple interpolation to — usually — a constant. Then by undoing the peeling-off 
transformations, one can determine very accurately the asymptotic expansion of the data 
up to a certain number of subdominant terms, which depends on the precision available 
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on the original data. This method has been applied to various nonlinear problems and 
shown to give very accura te expression of the a s ymptotic expansion when the data have 
enough precision (see e.g., IPauls fe Frischll2007t feardos et al]|2009l) . 
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